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Abstract 



We propose a form factor approach for the computation of the large distance 
asymptotic behavior of correlation functions in quantum critical (integrable) 
models. In the large distance regime we reduce the summation over all excited 
states to one over the particle/hole excitations lying on the Fermi surface in 
the thermodynamic limit. We compute these sums, over the so-called criti- 
cal form factors, exactly. Thus we obtain the leading large distance behavior 
of each oscillating harmonic of the correlation function asymptotic expansion, 
including the corresponding amplitudes. Our method is applicable to a wide 
variety of integrable models and yields precisely the results stemming from the 
Luttinger-liquid approach, the conformal field theory predictions and our pre- 
vious analysis of the correlation functions from their multiple- integral represen- 
tations. We argue that our scheme applies to a general class of non-integrable 
quantum critical models as well. 
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1 Introduction 



Form factors and correlation functions are central objects for the study of dynamical properties 
of models in quantum field theory and statistical mechanics. For general interacting systems, 
their calculation remains a fantastic challenge. In the integrable model situation in low dimen- 
sion, see e.g. [l- 14] and references therein, considerable progress has been made during the last 
forty years towards the exact calculation of these dynamical quantities 14I446J]. Of particular 
interest are the so-called critical models where the correlation functions of local operators are 
believed, from Luttinger liquid app roach and conformal field theory (CFT) predictions, to decay 
as power laws in the distance 47H57I]. There are several ways to approach the computation of 
correlation functions in such models, either from their massive field theories counterparts (in 
their short distance behavior) or from critical lattice models (looking at their long distance 
properties). In such limits one should recover the picture predicted by CFT, with the possi- 
ble additional insight into the explicit correspondence with the physical local operators of the 
microscopic initial theories. The latter information is crucial for obtaining full control of their 
dynamics. Indeed, while CFT is able to predict for example the algebraic coefficients in the 
operator product expansion of local operators together with the values of the various critical 
exponents driving the power law behavior of the correlation functions, it fails to determine the 
(non universal) corresponding amplitudes (vacuum expectation values and form factors) which 
are in fact deeply rooted in the detailed microscopic interactions of the physical models at hand. 

The purpose of the present article is to show that under quite general assumptions, which 
are valid in the integrable model situation, it is possible to determine exactly the asymptotic 
behavior of correlation functions in critical models starting from their expansion in terms of form 
factors. The use of form factors in this context might app ear to be a very strange strategy (see 
however 58J, |59(). Indeed, as has been demonstrated [6CH62I] . form factors of local operators in 



critical models scale to zero with a non trivial power law in the system size (or equivalently in the 
mass scale). This is not surprising as it just reflects the corresponding critical dimensions of the 
local operators (together with those of the creation operators of the excited states) considered. 
Hence the form factor expansion of correlation functions of local operators in critical models is 
a large (infinite) sum of small (vanishing) terms in the infinite volume limit. There it should 
be the collective (volume) effect of the sum that compensates the vanishing behavior of each 
term of this sum. Hence, extracting this collective effect directly from the form factor series 
seems a priori to be a quite involved task. However, a very interesting effect revealed itself in 
our recent derivation of the leading asymptotic behavior of correlation functions of the XX Z 
chain obtained from their re-summed multiple-integral representations [381 ] . Unexpectedly, it 
appeared there that the exactly computed amplitudes of the leading power law decay of the 
spin-spin correlation function were in fact given by properly renormalized (by a positive power of 
the system size) form factors of the local spin operator. The very appealing form of this answer 
was a strong motivation for us to develop a method for studying such correlation functions from 
their form factor expansion. 
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The method that we propose in the present paper opens up a way to analyze the large 
distance asymptotic behavior of correlation functions of critical models from their form factor 
expansion. Namely, it enables us to determine the leading behavior of each oscillating harmonic 
in the asymptotic expansion of the correlation functions, including not only the exact power 
law decays but also the associated amplitudes. To describe our method, we first study the 
correlation functions of quantum integrable models that are solvable by using the Bethe ansatz 
in their critical regime. The typical example considered in this paper is the integrable Heisenberg 
spin chain. However, the method that we describe in this paper is also applicable to the one- 
dimensional Bose gas (or quantum non- linear Schrodinger model), and more generally to any 
quantum inte grab le model possessing determinant representations for the form factors of local 



operators |33l. I63H65I]. In fact we will further argue that, under quite natural assumptions, it 
applies to a general class of non-integrable critical models as well. 

Generally speaking, independently on the integrability of a quantum model, the zero tem- 
perature limit of a correlation function of local operators located at two different space points 
x' and x' + x, say 0\{x') and 02 (x' + x), reduces to the ground state expectation value of their 
product: 

ww+x)) ^M^^^nm. (1J) 

Above, | ip g ) denotes the ground state and the coordinate x can be either a continuous parameter 
(for one dimensional field theories) or a discrete one (for lattice models). 

The form factor expansion of the zero temperature correlation function is obtained by in- 
serting a complete set of eigenstates | tp' ) of the Hamiltonian between the local operators: 



Thus the correlation function is represented as a sum of matrix elements (the so-called form 
factors) ^ ^(x) of the local operators Oj taken between the ground state | tp g } and any excited 
state \tp' ): 

What makes the representation (|1.2|) very appealing is its straightforward physical interpreta- 
tion. Moreover, it proved to provide a very efficient way to evaluate the correlation functions 



numerically as soon as exact representations for the form factors are known [461. I66H69J] . 

The form factor expansions are well understood in the case of massive models 1 914251] . 
Actually, in these models they are extremely effective in computing the long distance exponential 
decay of the associated correlation function. In particular, the leading long distance behavior 
of the correlation functions generally stems there from a few classes of excited states only. 
However, as we already discussed, the structure and properties of form factor expansions in 
massless (critical) models are more involved, as an infinite tower of classes of excited states 
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contributes to the leading asymptotic behavior of the correlation function, in contrast to the 
massive regime case. 

There has recently been important progress in computing the asymptotic behavior of correla- 
tion functions for quantum integrable models; these approaches rely basically on the Riemann- 
Hilbert analysis of related Fredholm determinants, and are quite technical 38l. l70l. |72|. |73(|. We 
believe that the method that we propose in this paper is simpler as it follows a path closer to 
the physical intuition and hence should be of a wider use (in particular beyond the integrable 
case) . 

As we just mentioned, for critical models, each term of the form factor series scales to zero 
as some non-trivial power 9 (depending on the excited state | ip 1 ) selected by the local operators 
considered) of the system size L [6fl (i'2 for L — > oo. In translation invariant systems we have 

• *P*.( x + x>) = L ~ e elxVex A (L4) 

Here the dependence on the distance of the form factor is contained in the phase factor e lxrpex , 
where V ex is the momentum of the excited state | ip' ) relative to the ground state. The amplitude 
A does not depend on x and remains finite in the thermodynamic limit. 

Our strategy for the analysis of the sum (|1.2[) is as follows. Initially, the sum (|1.2p runs 
through all the eigenstates and for a finite value of the distance x all eigenstates should indeed 
be considered. However, in the long distance regime, as follows from (|1.4|) . each term in the 
form factor series is quickly oscillating with the distance. This means that, similarly to what 
happens for oscillatory integrals, the leading behavior of the form factor series stems from 
the localization of the sum around specific critical points (like saddle points or edges of the 
summation interval). This fact allows us to carry out several reasonable simplifications of the 
summation range that do not affect the leading asymptotic behavior of the correlation function. 
First of all, we restrict the form factor sum (jl.2p to the summation over excitations of one 
specific type only, namely, the particle-hole excitations^. Second, in the large distance regime 
we take into account only contributions coming from excitations close to the endpoints of the 
Fermi zone, the so-called critical excited states 38|, [62j]. In this way, we naturally restrict 



ourselves to the excitations appearing in the conformal part of the spectrum of the model. We 
will show that the resulting restricted sum over critical excited states can be computed exactly. 
The asymptotic expansion for the two-point functions of the XXZ model that we obtain in 
this way confirms the predictions stemming from a correspondence with CFT/Luttinger liquid 
approaches and agrees with the part of the asymptotic behavior obtained through the exact 



Riemann-Hilbert methods 38l. I7CH73| 



The paper is organized as follows. In Section 2 we define the critical form factors and 
formulate our main strategy for evaluating their sums in the large distance limit. We give some 



1 Hence, in particular, bound states will not contribute in this approach, although we believe that they could 
appear in the analysis of dynamical correlation functions. Of course, for some models such as the non-linear 
Schrodinger model, all the excited states with finite excitation energy can be described in terms of particles and 
holes. 
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arguments confirming this approach. In section 3 we perform exactly the summation over the 
critical form factors. This computation is given in a general, model independent framework. It 
gives the leading asymptotic behavior of all the oscillating harmonics. In section 4 we apply 
this general result to compute the two-point functions for the Heisenberg XXZ chain. We 
reproduce the leading asymptotic behavior of the two-point functions as predicted by CFT 
and the correlation amplitudes obtained by the Riemann-Hilbert approach. Finally, in the 
conclusion, we argue that the above scheme applies to a large class of critical (integrable and 
non-integrable) models. In Appendix [A] we prove the summation formula which permits us to 
compute all the above critical sums exactly. 



2 Form factor series in the large distance limit 

Let us consider the case of quantum integrable models for which the ground state solution of 
the Bethe equations can be described in terms of real rapidities Xj densely filling (with a density 
/9(A)) the Fermi zone [— q, q}. In such models,the logarithmic Bethe ansatz equations for the 
ground state Bethe roots Xj take the form 0, y, 0-0] : 

N 

Lpo(Xj) ~ fl ^( A i " A *) = 2vr U - —5^ ) > j = l,-..,N. (2.1) 
k=i ^ ' 

Here L is the size of the model, and N is the number of quasi-particles in the ground state. 
The functions po(A) and $(A) are the bare momentum and bare scattering phase of the quasi- 
particles in the corresponding model. The value iV depends on the parameters of the model 
(such as the chemical potential or external magnetic field). When we take the thermodynamic 
limit L — > oo, the ratio N/L has a finite limit which we call the total density D = lirn^oo ^. 

In order to describe the form factor sum one needs to have a characterization of all the 
excited states as well. In the Bethe ansatz framework, the excited states are parameterized 
by rapidities {[ie a }¥' solving a set of equations similar to (|2.ip but involving other choices of 
integers l\ y . . . ,£j^i in the rhs: 

N ' ( N' + 1 A 
Lpo(Mj) ~ Yl " W k ) = 2tt ( lj — , j = l,...,N'. (2.2) 

k=i ^ ' 

Above, N' corresponds to the number of quasi-particles contained in the given excited state. 
The values of N' to be considered for carrying out the form factor sums depend on the type of 
local operator O that one deals with. A local operator O only connects states having a 'small' 
difference in the number of quasi-particles. In other words, for a given local operator, there 
exists an integer k such that only states with \N — N'\ ^ k give rise to non-zero form factors. 
For instance, later on, we will focus on the XXZ spin-1/2 chain. There, we will consider the 
local spin operators, where either N' = N or N' = N ± 1. We do however stress that our 
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method allows one to consider more complicated cases of local operators that connect more 
distant values of N and N' . 

For excited states of the type that we are interested in, it is convenient to describe the 
solutions of the Bethe equations in terms of a set of integers alternative to £i, . . . ,£n'- Namely 
one characterizes the excitations in terms of holes in the Dirac sea of integers for the ground 
state £j = j and additional, particle-like integers outside of the zone {1, . . . , N'}. In other 
words, the integers lj describing an excited n-particle-hole state take the form 

l a = a , ae{l,...,N'}\{h 1 ,...,h n } £ ha = p a , p a € Z\{1, . . . ,N'}. (2.3) 

Each choice of pairwise distinct quantum numbers {p a } and {h a } characterizes an excited state 
associated with a configuration of the particle and the hole rapidities {^ Pa } and {nh a }, see e.g., 

The particle/hole excitation contribution to the two-point function has the form 

(O x {x') 2 (x + x')) ph = lim L ' 6 ^ A (M> KIIW, {h})- (2.4) 

{P},W 

Here we have explicitly indicated that the amplitude A (|1.4|) depends on the particle/hole 
rapidities {// p } and {/x^} as well as on the corresponding integer quantum numbers {p} and {h}. 



Let us comment briefly on this dual dependence (see 62] for more details). If all particle/hole 
rapidities are separated from the Fermi boundaries, then the amplitude smoothly depends only 
on these quantities A = A({fj, p }, {nh})- However, as soon as these rapidities approach the Fermi 
boundaries, the amplitude shows a discrete structure: a microscopic (of order 1/L) deviation of 
a particle (hole) rapidity, corresponding to a finite shift of the integer quantum numbers {p} or 
{h}, leads to a macroscopic change of A. Hence, in this case, one should take into account the 
discrete structure of the amplitude, and therefore we write A = A({fi p }, {[J>h}\{p}i {h})- 

In the following, we will argue that the oscillatory character of the sum in (|2.4|) localizes 
it, in the absence of any saddle points of the oscillating exponent^, in the vicinity of the Fermi 
boundaries ±q. For this purpose, we introduce several definitions [62|. 

Definition 2.1. An n particle-hole excited state {^e a } is called a critical excited state if the 
rapidities {/x Pa } and {nh a } defined by such a state accumulate on the two endpoints of the 
Fermi zone in the thermodynamic limit. Form factors corresponding to any such a state are 
called critical form factors. 

The critical excited states are characterized by the specific distribution of particles and holes 
on the Fermi boundaries, ±q. Namely, assume that, in the thermodynamic limit, there are 
particles whose rapidities are equal to ±q and holes whose rapidities are equal to ±q. Then, 
obviously, 

n t + n v =n h+ n h = n - ( 2 - 5 ) 



2 Such saddle points could appear in the time dependent case and will be considered in a forthcoming publi- 



cation. 
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Definition 2.2. A given critical excited state belongs to the class if the distribution of 
particles and holes on the Fermi boundaries is such that 

n p~ n h =n h~ n p = ~ n ^ ^ ^ n - (2-6) 

Then such an excited state has momentum 2£k F in the thermodynamic limit, where k F = ttD 
is the Fermi momentum. Its associated critical form factor will also be said to belong to the 
class. 

Now we would like to argue in favor of the localization of the form factor sum on the Fermi 
boundaries in the large distance limit. For this purpose we will use an analogy with oscillatory 
integrals. Hence, suppose that we deal with some multiple integral of the type 

in(x)= J j ^/(w.win^ 1 ^ 1 ' X ^°°> ( 2 - 7 ) 

R\[-g,«] [~q,q] j=1 

where /({/%>}> {t^h}) is a holomorphic function in a neighborhood of the real axis and p(n) has 
no saddle points on the integration contours. By deforming the integration contours into the 
complex plane, one makes the integrand exponentially small everywhere except in the vicinity 
of the endpoints ±q. Then the large x asymptotic analysis of (|2.7p reduces to the calculation 
of the integral in small vicinities of the endpoints, where /({/%,}, can be replaced by 

f({±q}, {±q}). Indeed, such replacement does not alter the leading asymptotic behavior of 
I n (x). Now, if f({fj, p }, {/J-h}) has integrable singularities at ±q, for example /({yU p }, {^h}) = 
(q — fih 1 Y + (fJ'h 1 + q) v ~ freg({fJ> P } , {fJ>h})i then, when carrying out the asymptotic analysis in the 
vicinities of ±q, one has to keep the singular factors (q =F Hhi) u± as they are, but we can replace 
the regular part freg{{fJ> p }, {fJ>h}) by an appropriate constant f reg ({±q}, {±q}). Again, this 
approximation does not alter the leading asymptotic behavior of I n (x). 

Thus, working by analogy, we may fairly expect that, as soon as V ex has no saddle point, the 
sum (|2.4|) will localize on the endpoints of the summation for the integers p a and h a , namely, 
in the language of rapidities, on the two endpoints of the Fermi zone. In other words, the sum 
will localize on the critical excited statefi The part of the amplitude A({fj, p }, {fJ>h}\{p}i {h}) 
smoothly depending on the rapidities plays the role of the regular part f reg in the multiple 
integrals (|2.7p . There one can set {n p } and {^h} equal to their values in the given P^ class: 

MM, M\{p}, {h}) — ► A({q} + U {-q} {q}, U {-q} \{p}, {h}) = A^({p}, {h}). 

V P h h 

(2.8) 

Such replacement should not change the leading asymptotic behavior over x of the form factor 
series. 

On the other hand the remaining part of the amplitude A^\{p}, {h}) explicitly depending 
on £ and on the quantum numbers p and h plays the role of the singular factors (q =p ^h 1 ) v± 

3 Note that for discrete sums, the integration by part procedure could be carried out, leading in simple cases 
to such a result in a perfectly controlled way. 
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in the integral fj2.7f) since, for large system size L, it varies quickly in the vicinity of the Fermi 
boundaries. In fact, in the thermodynamic limit, this rapid variation of the amplitude as a 
function of the rapidities manifests itself as a kinematical pole of the form factor whenever one 
particle and one hole are approaching simultaneously the same Fermi boundary. In terms of 
the quantum numbers {p} and {h} however, there is no singularity and the coefficients Ay) 
are always well defined. So, fixing the class of critical form factors, we should take the 
sum involving A^({p}, {h}) over all the excited states (namely over all the possible quantum 
numbers {p} and {h}) within this class before taking the thermodynamic limit. 

This picture fits perfectly into the scheme of approximations used in the application of CFT 
(or Luttinger liquid) methods to the computation of the asymptotic behavior of correlation 
functions. Moreover, it gives a clear interpretation of why these predictions actually work and 
are so universal. 

Let us point here an important remark. Although we just used the multiple-integral asymp- 
totic behavior analogy, the sum that we consider here cannot be replaced straightforwardly by 
an integral sum even for large system size L. The non-integer nature of the critical exponent 9 
would make such a correspondence quite subtle, eventually producing divergent integrals times 
remaining vanis hing factors of L. This feature explains the difficulties already noted in the lit- 
erature, see e.g. 744761] . while trying to use the form factor approach directly in the continuum 
limit for massless models. The main advantage of our approach (in comparison to a field theory 
framework) is in computing exactly the form factor sum for L large but finite, hence producing 
an explicit compensation for the vanishing factor L~ e and making the thermodynamic limit 
L — > oo straightforward, as we will show in section O 



3 Summation of critical form factors 

We have already insisted that critical states are gathered into various classes having distinct 
values of their excitation momentum in the thermodynamic limit. More precisely, all states 
belonging to the class have momentum 2lk F . Also, we stress that all critical states have, in 
the thermodynamic limit, a vanishing excitation energy (the latter goes to zero as 1/L), this 
independently of their class. 

Actually, similarly to what happens for the excitation energy, for L large but finite, the 
momentum of any critical state deviates from 2lk F by terms of order 1/L. Namely, it follows 
from the Bethe equations (|2.1|) and ([2.2|) that 

N' N 2ir n 

Vex = ^PoiVij) ~ ^Poi^j) = -JJ^2(Pk ~ h k), (3-1) 
j=l j=l k=l 

where p and h are the integer quantum numbers of particles and holes. It is convenient to 
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re-parameterize these integers for the critical excited states as follows 

Pj =Pj + N\ if n Pj = q , 

hj=N' + l-h+, if ii hj =q, { ' 1 

hj = hj, if fj, hj =-q. 

All the integers {p^} and {h^} defined above are positive and vary in a range such that 

lim i-Vp±= lim ly^± = 0. (3.3) 

It means that all the particles and holes collapse to ±q in the thermodynamic limit. Then the 
expression for the excitation momentum of critical states belonging to the class becomes 



V ex = ^W' + ^V«\ (3.4) 



Above, we have set 



j=i j=i j=\ j=i 

We recall that in the thermodynamic limit N'/L — » with D = k F /TT. 

It is customary in the Bethe ansatz framework to describe excited states in terms of their 



shift function, see e.g. [12| . The latter characterizes the way in which the Bethe roots (rapidities 
of quasi-particles) of a given excited state are shifted with respect to the ground state ones. It 
is in fact the result of the non trivial interactions that hold in the model, making the quasi- 
particles inside the Fermi zone move slightly due to the creation of particle/hole excitations. In 
the case of a generic particle/hole excited state, this shift function depends on the rapidities of 
the particles {fi p } and holes {/J-h} occurring in the given state. However, in the case of critical 
excited states, this shift function is the same for all the states belonging to the same class, 
and is denoted by Fi(X). Of course, the function Fe(\) depends on the model, and on the 
excited state selected by the local operator under consideration. For integrable models, it is 
obtained as a solution of a linear integral equatiorjf). It was shown in 6C , 62 ] that the exponents 
9 can be expressed in terms of the Fermi boundary values of this function: 

Ff = F t (-q), F+ = F e (q) + N' - N. 

In particular the exponents 6 for the form factors of the class P^ are all equal to a same value, 
6i, which can be written as [62| 

9 e = (F- +£) 2 + (F+ + £) 2 . (3.6) 



4 The linear character of the integral equation satisfied by the shift function implies that Fe(X) linearly depends 
on the integer I (see section Q . 
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As was mentioned above the values of the shift function F» depend on the specific model and 
on the operators Oj that we deal with. However the functional expression (|3.6p for Og in terms 
of the constants is universal and model independent. 

Thus, the sum over form factors (|2.4p being restricted to the critical ones takes the form 



(0 1 (x')0 2 (x + x')) cr = lim y L- e <e 2ixek r V e^'iW^},^}). (3.7) 

l= — O0 {P ± },{h ± } 



Observe that this sum is separated into two parts: the external sum over different classes, 
and the sums over the quantum numbers {p ± }, {h^} within a fixed class. We call the last 
ones critical sums of order t. Calculating these critical sums we deal only with the discrete 
amplitude AM^},^}) weighted with the exponents exp(^p7e:c )• If is remarkable that 
A" has a purely kinematical interpretation and that, up to a constant factor, its functional 
form is universal and model independent. 

Consider the ^-shifted state | tpi ) with the Bethe roots given by the following equations: 

N ' f N' 4- 1 \ 

LPo(j*j)-J2#(H-l*k)=2'x[j+t j = l,...,N'. (3.8) 

k=l ^ ' 

This excited state belongs to the P^ class. Since all the form factors of the Pi class scale as 
L~" e , we define the renormalized amplitude F^F^ which corresponds to matrix elements of 
the operators 0\ and O2 between the ground state and the ^-shifted state | ): 

T (1) F (2)_ ^ (^101(^)1^X^102(^)1^) ,, q . 



This quantity is finite in the thermodynamic limit 611,162] and, due to the translation invariance, 
it does not depend on x'. 

The discrete amplitudes corresponding to a given P^ class can be expressed in terms of the 
above quantity. Namely, one has [62] 



A^{{ P ±},{h±})=F^F^ 



G 2 (l + F+)G 2 (1 - F7) /sin(7rF+) V"" /sin(7rF, 



G 2 (l + e + F+)G 2 {l-£-F-) V tt J \ n 
x Rni,ni({P + }dh + }\F+) R--({p-},{h-}\-F e -), (3.10) 



where G(z) is the Barnes function satisfying G(z + 1) = T(z)G(z), and 

n n' 

u( P] - p k ? nth, - h k f n n , 

R n y({p}, {h}\F) = ——^ U U . (3.11) 

n u(pj+h k -i) 2 fc=i vw k = i { ' 
3=1 k =i 
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Let us comment the formulas above. 

The representation (|3.10p was derived straightforwardly for the particle-hole form factors in 



the XX Z Heisenberg chain [62] and quantum one-dimensional Bose gas 721 ] . It represents the 
thermodynamic limit of known determinant formulas for form factors for finite N and L 33j, [771] . 
In this representation, only the constant J-^J 7 ^ depends on the model and on the operators 
considered. The rest of the equation (|3.10p is universal. This occurs because the particle-hole 
form factors for both models are proportional to generalized Cauchy determinants, namely, 

N N' 

n( A i ~ x k) milt, -/%) 

JTIP • ^ 3 - 12 ^ 

n n(Ai-^ fc ) 

3=1 k=l 

We call the rhs of ()3. 12|) the generalized Cauchy determinant because for N' = N it gives known 
expression for the standard Cauchy determinant det[l/(Aj — (J>i k )] in terms of the rapidities \j 
of the quasi-particles in the ground state and fj,£ k in the excited state. The appearance of such 
an object in the formulas for form factors has a clear physical interpretation. The fact that 
•Fipgi/j' = as soon as Xj = or fj,g. = /j,g k simply takes into account the fermionic structure 
of the particle/hole description of the ground state and of the excitations. The singularities 
at Xj = fj,£ k produce, in the thermodynamic limit, poles occurring whenever the rapidities of 
particles and holes coincide at the Fermi boundaries. Such a behavior near the Fermi boundary 
has a purely kinematical origin and is quite universal. It appears for example in the form factor 
bootstrap program. There the residues of the n-particle form factor at the particle/antiparticle 
poles are given in terms of (n — l)-particle form factors [2CJ 



On the other hand it is not difficult to check (see e.g. [62J) that the combination of the 
generalized Cauchy determinants corresponding to the ratios (^l^/^i/w, of form 

factors gives, in the thermodynamic limit, exactly the universal part of the amplitude (|3.10|) . 
provided J 7 ^^, and J~^\ are critical form factors of the class. Thus, a rather complicated, 
at first sight, dependence of the form factors on the quantum numbers such as in (|3.10p has a 
simple physical origin. 

One can observe that the expression (|3.1U|) is factorized into the product of two parts corre- 
sponding to the left and right Fermi boundaries. Such factorization occurs in the thermodynamic 
limit, when we neglect all corrections of order o(l) at L — > oo. If we do not neglect these correc- 
tions, then we evidently produce subleading contributions to the form factors of order L — <— 1 , 
£-0<-2 e ^ c Q n th er hand in this case we do not have a complete factorization of the final 
answer, like in (|3.10p . and thus, we can say that higher order corrections describe an effective 
coupling between the two Fermi boundaries. 

Using the representation (|3.10p we can present the sum of the critical form factors in the 
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following form: 



(0 1 (x')0 2 (x + x')) cr = lim V L- e '^ k *jf)jf)-— 

D 

2irix 



G\l + F;)G\l-F t 



+ t + F+)G 2 (l-£-F; 



£ exp 

{ P ±},{h±}^ 

P h h P 



L 



sin(7rF/)\ 2n " /sin(^~)\ 2nh 



7T 



7T 



X ^n+({P + H^ + }l^ + ) ^-„^({P-}.{^"}I -^")- ( 3 - 13 ) 

We would like to mention that this equation gives explicitly the contribution of the critical 
form factors in the thermodynamic limit. 

It is important to notice that the critical sums like in (|3.13p are universal and hence appear 
already in the models equivalent to free fermions. There, of course, the shift function is a trivial 
constant. Indeed, in free models one cannot obtain a non-trivial shift function because of the 
absence of interaction. Therefore in such models a shift between the ground state quasi-particles 
and the ones describing the excited state either occurs due to non- vanishing difference N — N', 
or it can be created artificially by considering twisted excited states (see section In contrast, 
for the interacting models we deal with, a non-trivial shift function arises automatically. It is 
remarkable, however, that being restricted to the states belonging to a fixed class, this non- 
trivial shift function enters the critical sums only through its values at the Fermi boundaries, 
namely through the two constants F» . In this sense one can say that we reduced our model to an 
effective (deformed) free fermion one in every class. We stress however, that for interacting 
systems such a reduction cannot be carried out uniformly for all critical form factors, in contrast 
to what happens for a free theory: every P^ class of form factors should be described by its own 
(deformed by the values Ft ) free fermion theory. 

The most important property of the critical sums is that they can be explicitly computed 
due to the following summation formula: 



E E E 

n.n'^O Vi< — <Pn hl<—<h„ 

n-n'=i p a eN* h a £N* 



exp 



2nix 
L 



5> 

i=i 



1) + 



k=l 



sinirF 



G 2 {l + £ + F) 
G 2 (l + F) 



7T 



exp 



2n' 



Rn,n>({p},{h}\F) 



exp j 



(3.14) 



We prove (j3.14j) in Appendix A. Amazingly, for the case I = 0, such a summation formula 
appeared already in the problem of defining Z-measures on partitions 78-84 f|. Note that this 
formula was also useful in the case of thermal correlation functions of the Bose gas 85 ] . 



5 We are grateful to O. Lisovyy for pointing out these works to us. 
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Using this identity we obtain for the two-point function 

{Oi(x')0 2 (x + x')) cr = lim V t e - l - . (3.15) 

£=-oo f 1 _ e — ) f 1 _ e ~ — 



Now the thermodynamic limit can be easily taken using (|3.6p and we obtain our final result: 

(0 1 (x')0 2 (x + X% r = y e ^k F+l Zsi g n( X ){(F- + l)*-(F+ + e)l] S e J\ (316) 

^co (2vr|x|)^ 

Observe that the critical exponents 0£ driving the leading power law decay in (|3,16p coincide 
with the ones describing the scaling behavior of the form factors (in terms of the system size). 
This property is quite natural as the distance x and the system size L in (|3.14p appear only in 
the combination x/L. As soon as we expect the sum over critical form factors to compensate 
the vanishing coefficient L~ 9e , it should produce a factor (x/L) £ , leading to an answer like 
in (|3,16p . The 1/L-corrections to the form factors describing the coupling of the two Fermi 
boundaries, and mentioned above, should produce respective 1/x-corrections to the leading 
asymptotic behavior of each oscillating harmonic in (|3. 16[) . In the CFT language such higher- 
order corrections should correspond to the contribution of the descendants of primary fields. 

We would like to underline once again that this is an explicit result for the contribution of 
the critical form factors to the two-point functions of a very large class of integrable critical 
models. Our main conjecture states that this equation gives the leading asymptotic behavior 
of all the oscillation harmonics of the two-point function. This conjecture is supported by the 
fact that, in the case of the XXZ chain (or the quantum non-linear Schrodinger model), we 
reproduce the exact answer stemming from a Riemann-Hilbert problem analysis for the first 
few harmonics. Also, our asymptotic formula reproduces the CFT predictions for the XXZ 
spin chain and quantum non-linear Schrodinger model for all the oscillating harmonics. 



Our method gives a simple explanation of the main result of 38|, [6l( . We have shown there 
that the amplitude of the leading oscillating term of the spin-spin correlation function for the 
XXZ spin chain (or of the density-density correlation function for the quantum non-linear 
Schrodinger model) can be expressed in terms of a rescaled special form factor. Moreover we 
observed that the critical exponent driving this leading oscillating term is equal to the one 
used for the rescaling of this special form factor. The critical form factor approach gives a 
simple explanation for both results, generalizing it to the leading asymptotic behavior of all the 
oscillating harmonics. Recently another explanation of these facts was given in (8f|. 

Let us conclude this section by one remark about the non-oscillating (£ = 0) term in (|3.16j) . 
For some correlation functions the corresponding exponent is zero 9q = 0; thus the corresponding 
term is constant. A small modification of the critical form factor approach however permits 
us to compute the first non-oscillating correction to this formula. It requires the notion of a 
'twisted propagator' and 'twisted Bethe states'. In section0]we will show how these objects can 
be introduced for the XXZ spin chain. 
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4 Correlation functions of spins in the XXZ chain 



In this section, we apply the general method described above to the spin-spin correlation func- 
tions of the XXZ spin-1/2 chain. Recall that this model is given by the Hamiltonian 

L , L 

# = £ KX+i + <r V A +1 + A(o{o£u - 1)) - ^ E < C 4 - 1 ) 

k=l k=l 

Here cr^.' y ' z are the spin operators (Pauli matrices) acting on the k-th. site of the chain, h is an 
external magnetic field and the model is subject to periodic boundary conditions. The XXZ 
spin chain above exhibits different phases depending on the value of the anisotropy parameter 
A. Here, we only focus on the critical regime — 1 < A < 1 where we set A = cos£. 

In the following, we consider the two correlation functions^ (o"fo"m+i} an< ^ ( a i a m+i)i Wi th 
ct ± = \{o x ± ia y ). Note that, since the lattice spacing is discrete, we have denoted it by m 
(instead of x as in the previous sections). 



4.1 The correlation function (cricr m+1 ) 

The correlation function (o"J f <T~ +1 ) can be computed using the general method. 

The only excited states which produce non-zero contributions here are the Bethe states with 
N' = N + 1. The critical values of the shift function were computed in 62]: 

Ff=£(Z-l)-± F+=£(Z-l) + ± (4.2) 

where Z = Z(±q) is the Fermi boundary value of the dressed charge Z(X) given by the following 
integral equation: 

q 

Z(X) + — fdfi — — S ' n2C , ZU) = 1. (4.3) 

v ; 2vr J p sinh(A - fi + i() sinh(A -fi-i() 



The critical exponents are given by 



6 e = 2e 2 Z 2 + ^. (4.4) 



We introduce the simplest form factor of the class: 

The rapidities corresponding to the ^-shifted state solve the following logarithmic Bethe equa- 
tions: 

N+1 ( N\ 

Lpo(ti j )-Y,^H-^k) = 27T(j + £ + l--\, j = 1, . . . ,N + 1. (4.6) 

k=i ^ ' 



'{a 1 cr+ +1 ) can be obtained from {at^m+i) changing m into 
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The bare momentum and the bare scattering phase for the XX Z model are 



p (A) = «lo£ 



'sinh(f + A)* 
sinh(f - A) 



tf(A) = ilog 



sinh(i£ + A) 



v sinh(i£ — A) 

Applying the general formula (|3.16p . we obtain the asymptotic result 

( i \ m — 1 



2™^ 



l^limlk 



* \T1 



(27rm) 



2£ 2 .Z 2 • 



(4.7) 



(4.8) 



The dominant term of this expansion (t = 0) behaves as predicted by the Luttinger liquid 
approach and CFT. 



4.2 The correlation function (ofcr^ +1 ) 

In order to compute the two-point function (o"fcr^ +1 ), it is convenient to introduce its generating 
function ( e 2 ™ 2 ™) 0, 0, Q , where 



(4.9) 



k=l 



This function generates the two-point function (crfa^ , x ) via 



o=0 



21? + 1, 



(4.10) 



where is the second lattice derivative. 

The large m behavior of {e 27UQ ^ m } can be obtained in the framework of the form factor 
approach. For this one introduces the complete set of eigenstates of the twisted transfer matrix 



T K {X) (see e.g. [371] ). where k = e is a complex (twist) parameter. This transfer matrix 
produces an XX Z spin-1/2 Hamiltonian but subject to quasi-periodic (instead of periodic) 
boundary conditions. Hereby the excited states of the Hamiltonian (I4.1ja re obtained as the 

881 ] . The introduction 



37, 61 



ol — y limit of the eigenstates of the twisted transfer matrix T K (A) 
of the extra parameter k is very useful in the analysis of the form factor series. 

In complete analogy with the standard algebraic Bethe ansatz considerations, the eigenstates 
\^k{{iA)) °f T~ K (\) can be parameterized by sets of solutions of the twisted Bethe equations 

N 



£po(Mi)-5^0(Mj -Mfc) 



k=l 



2-7T ( rij + a 



N + 1 



J 



1, 



,N. 



(4.11) 



For this correlation function only the form factors with N' = N produce non-zero contributions. 



7 There is an extra general ( — 1)'" factor, which can be removed by the re-definition of the Hamiltonian. 
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From the solution of the quantum inverse problem it can be easily shown that for any 
eigenstate of the twisted transfer matrix the form factors are proportional to the scalar products 



<<MM) \e 2maQ -\^ g ) = e^™<^({//}) (4.12) 

where V ex is the momentum of the excited state. 

Critical form factors related to the operator e 2nia ^ m were computed in 62]. For the 
class, the boundary values of shift function Fi are 

Ff = F+ = aZ + l{Z-l). (4.13) 

Thus the exponent 8e(a) can be written as 

e (a) = 2^(a + £)Zy. (4.14) 
The amplitudes in this case are just the normalized scalar products 

\T J 2 = lim L 9 'M . , l^f'tHlC m ( 4 - 15 ) 
where | ip' a+ i ) is the state parameterized by the Bethe roots of 

Lp (»j)-J2' & (»3-Vk) = ^(j + a + e , j = l,...,N. (4.16) 

k=i ^ ' 

Note that ( ip g \^' a+ () = at a = and I ^ 0, as the scalar product of two different eigenstates. 

Now, using the general result (|3,16p . we obtain the leading asymptotic terms for all the 
oscillating harmonics: 

\ lor | a+t \ (2vrm)^( a ) 



l=-oo 



The dominant terms with I = 0, ±1 were already obtained by the Riemann-Hilbert analysis 
of the multiple-integral representations 38|, [6l|. Here, we get in a much simpler way the same 
result as well as the natural extension to other harmonics reconstructing the a-periodicity of 
the correlation function. 

To obtain an asymptotic expansion for the two-point function we should take the a-deriva- 
tives and lattice derivatives of (|4,17p . Here one should remember that the coefficients l-T^+J 2 
with I ^ have a second-order zero at a = 0, while the coefficient | T a | 2 does not vanish in this 
point. This leads to the following result: 

2Z 2 °° 1 
(<J z a z m+1 ) cr = (2D-l) 2 --^ + 2^cos(2m£fc F )|^| 2 - ^ 2 , (4.18) 



(2irmY 
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where 



lim L 

L— KX) 



2£ 2 Z 2 



\{i> 9 \° z M)f 



(4.19) 



and the Bethe roots corresponding to the ^-shifted state are defined by (|4,16p with a = 0. 

The first terms of this expansion (non-oscillating terms and 1 = 1 term) reproduce the 
asymptotic behavior predicted by various methods (Luttinger liquid approach 47H50I]. CFT 
5ll455j|). The amplitude of the first oscilla ting term is given by the special form factor as 



predicted by the Riemann-Hilbert approach [38|, [61] . 

We would like to underline that from the beginning we could apply the general method 



directly to the two-point function without considering the generating function (e 



2-iriaQn 



). How- 



ever the method proposed here (the use of twisted eigenstates) has one evident advantage: it 
produces the first subdominant non-oscillating term of order —j (the dominant non-oscillating 
term being a constant). 



Conclusion 

In this paper we have developed a new method for computing the contribution of the critical 
form factors to the asymptotic behavior of two-point correlation functions of critical integrable 
models. This result is quite general and can be applied to a very large class of integrable field 
theories and lattice models. We have shown that for the Heisenberg spin chain this contribu- 
tion reproduces the CFT predictions for the asymptotic behavior of the correlation functions. 
Exactly the same result can be obtained for the quantum non-linear Schrodinger model, repro- 
ducing once again existing predictions for the asymptotic behavior. 

We hope that the method that we presented here can be applied to a wide class of non- 
integrable systems as well. We have pointed out already that, within every fixed class, the 
critical sum over form factors coincides with the one corresponding to an effective (deformed) 
free fermion model. Therefore, if we have a general quantum system which admits, for every 
fixed excitation momentum, a (deformed) free fermion interpretation of its low energy spectrum, 
we can expect the sum over form factors within this sector of excitations to coincide with the 
critical sums considered in this paper. Then, calculating these sums explicitly via ()3.14j) . one 
obtains the leading critical exponents for all oscillating harmonics in the asymptotic expansion 
of the correlation function. 

We would like to stress once more that effective free models describing the critical form 
factors of interacting systems are different for different classes. In particular, the values of 
the shift function on the Fermi boundaries change from one class to another. Nevertheless, 
these values are certain constants for all excited states of the same class. It is quite natural to 
expect that, in the thermodynamic limit, these constants will only depend on the P^ class of the 
excited states, and not on other details, in much the same way as in the integrable situation. 
In fact, the shift function just encodes the polarization of the Dirac sea right at the Fermi 
boundaries as the result of the microscopic interactions at work in the system considered. Then 
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all the machinery that we described for the integrable situation will lead to the computation of 
the asymptotic behavior of the correlation function. The only difference comes from the fact 
that, in that case, the shift function F does not satisfy some exact linear integral equation and 
that its value at the Fermi boundaries could be hard to compute. However this quantity should 
be just considered as a physical parameter to be plugged into the summation formula ()3.14j) 
to determine the corresponding asymptotic behavior of the correlation function. This picture 
gives a nice (fully interacting) explanation of the origin of the critical exponents: they are just 
the result of the polarization of the Dirac sea at the Fermi boundaries produced by the creation 
of particle-hole excitations. Moreover, there is here some notion of universality. Generically 
different models possess different shift functions. However, models having the same values Ff~ 
will have the same critical behavior. We believe that, all together, this scheme provides a 
microscopic particle-type description of CFT starting from lattice models that should be useful 
both from conceptual and practical viewpoints. 

In a forthcoming publication we will show how this method can be extended to study the 
time dependent correlation functions. 

Acknowledgements 

J.M.M., N.S. and V.T. are supported by CNRS. N.K, J.M.M. and V.T. are supported by ANR 
grant ANR-10-BLAN-0120-04-DIADEMS. N.K. is supported by the Burgundy region, FABER 
grant 2010-9201AAO047S00753. We also acknowledge the support from the GDRI-471 of CNRS 
"French-Russian network in Theoretical and Mathematical Physics" and RFBR-CNRS-09-01- 
93106L-a. N.S. is also supported by the Program of RAS Mathematical Methods of the Non- 
linear Dynamics, RFBR- 11-01-00440, SS-8265.2010.1. When this work was done, K.K.K. was 
supported by the EU Marie-Curie Excellence Grant MEXT-CT-2006-042695 and DESY. N.K., 
N.S. and K.K.K. would like to thank the Theoretical Physics group of the Laboratory of Physics 
at ENS Lyon for hospitality, which makes this collaboration possible. N.K., J.M.M. and V.T. 
would like to thank LPTHE (Paris VI University) for hospitality. 

A Proof of the summation formula 

In this appendix we give the proof of the formula for the critical sums (|3.14p . We start with the 
simplest and most fundamental case £ = 0. Given < 1 and a complex valued is, we define 
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the function fo(v,w) by the series 



ac.«) = E E E ^.^-"(^) 

x ( det „^^) 2 n ^fe|^(;,-^) , (A „ 

V Pj + h k -ij ax r 2 (p fc )r 2 (/i fe ) 



n=0 pi<— <Pn hi<---<h n 

p a ew h a eN* 



The summand is a symmetric function of the variables {p} and likewise a symmetric function 
of the variables {/i}. Moreover it is vanishing whenever two of the p-type variables or two of 
the /i-type variables coincide. Therefore, 



E E =dy> E E ■ («) 

pi<— <p« hx<---<h n pi,...,p n =\ hl,...,h„=l 

PaGN* ft. aG N* 



Notice that, given any symmetric function {/t}) of the {p} and {/i} variables, one has 

oo oo / \ 2 

E E (^j-tt^t) 

pi,...,p„=l hx,...,h n =l v 7 
This property allows one to recast the series (jA.ip as 

OO j oo 

/ (z,,u;) = V — V det^i nV(hj,h k ), (A.4) 

' * 77-1 fc=l,...,n 



n! 

n=0 ftx,. ..,/ln=l 



where 

™ ^_v^ ^ +(fefc+ ^ )/2 r 2 (p + i + ^) VsinvT^V r(fe fc -^)r(^-^) 
^' * J "^V(p + fc*)(p + ^) ^(p + i) ){ vr j r(fc fc )i%) • ( 5) 

By definition, (jA.4j) gives the determinant of the following infinite matrix: 

f (u, w) = dct , : « [5 jk + V(J, k)]. (A.6) 

k— 1, . . . ,oo 

This determinant can be computed for general v and w. In fact, whenever v is integer, this 
infinite determinant boils down to a finite size determinant which can be computed using the 
Cauchy determinant formula. In the general case, however, no such simplification arises and 
one has to build the computation on known determinant evaluations of Toeplitz and Hankel 
operators obtained by H. Widom 89]. Below we consider the two cases separately. 
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A.l Integer v 

Making the replacement hj pj, we conclude that fo(u,w) = fo(—v,w). Therefore it is 
sufficient to consider the case v = N, where N € Z+. 

If v = N, then the combination sin 2 (7ri/) T(hk — v)Y(hj — v) vanishes as soon as hj > N or 
hp- > N. Thus, the series in each hj in (|A.4|) becomes truncated at hj = N, and we obtain the 
determinant of an N x N matrix: 

fo(u,w) = dot ... : x [6 jk + V(j, fc)] . (A.7) 

k=l,...,N 

Setting v = N in (|A.5[) we get, after simple algebra, 

TV 

u; (j+fc)/2 °° n(p+^) 2 



Wt) = nS V W P " =1 , W r. (A.8) 

1 1 u ~~ m ) 1 1 (A ~~ m ) p 

m— 1 m — 1 

Let us transform the det(J + V) obtained as 

, lr „. det N ( AA T + AVA T ) .. . 

with 

Aj k = w~ k l 2 W-\ j,k = l,...,N. (A.10) 

Obviously, 

N N-l 

(det at A) 2 = W ~ N ( N +^ 2 H(j - kf = w~^ N +m H (k\) 2 . (A.ll) 

j>k k=0 

To compute the product AVA T one can use the following identity: 

£ -R = (-1)^^ > , = l,..,iV, (A.12) 

e=1 (£ + P) IT it - m) U(P + m) 

m=l m=l 

which follows from 

/z s ~ l dz 
= 0, s = l,...,N, with R > max(A^,p). (A. 13) 

|2|=jj (z + p) l\ (z- m) 

m=l 



Then we obtain 



w N(N+l)/2 r N 



fo(v,w) = -jj-^ det at 

n (kir 

k=0 



w- p pi +k - 2 + (-iy +k ^V +fe ~ 2 

p=l p=0 



(A.14) 
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Setting here w = e , we get 

e -tN(N+l)/2 



N-l 

n (*o a 

fc=0 



det 



N 



N 



4 



+k-2 



-pt 



p=i 



p=0 
~ N-l 

n w 

fc=0 



det 



JV 



+Ai-2 e 



JVt 



1-e- 



(A.15) 



We now use that 



lim 



detjv*(nj,u fc ) t-/ 1 

= llmW det ^ 



"li-..i«jy->«o 



fc=0 



(A; 



0* 



+fc-2 



$(u ,uo) 



(A.16) 



where u) is any two-variable function such that all derivatives entering (|A. 16|) exist. The 
notation A means the Vandermonde determinants of the corresponding variables. 

Clearly the equation (|A.15p can be written as a homogeneous limit of the type (|A.16D : 

r „N(t+uj+v k ) 
M»,e-*) -' V ' V -'"-' ^ ' ' ' 



-tN(N+l)/2 Um det 



N 



,...,v N -iO 



1 _ e -t-Uj-v k 



£sr\u)£sr\v). 



(A.17) 



It is a Cauchy determinant; hence 



lim ^gn^ v 



ui,...,« w ->o A(u)A(t>) 

«1 DM-tO 



nc 

j,fc=l 



-t-Uj-V k 



r 1 - 



(A.18) 



(A.19) 



The calculation of the homogeneous limit becomes now trivial and we obtain 

fo(u,w) = (l-w)- N \ 

A. 2 Non-integer v 

Here we only give a sketch of the proof. Let us have \w\ = 1 but w ^ 1. Let also v satisfy 
— 1/2 < 5R(z/) < 1/2. The general value of v can be reached by analytic continuation. 

It is readily seen that the infinite matrix V in (|A.6p can be factorized as a product of two 
matrices: 

f (v,w) =det [I-U{p)U{-p)\, 

where 

,(j+*-i)/a r(j - v)Y{k + i/) 



Uj k {v) 



simru 



n{j + k — i) r(j)r(fc) 

Let and T(i/) be the following Toeplitz matrices: 

w (h-j)/2 



j — k + v 



and 



j,k = l,2,.... 



T jk (v) = wi- k T jk (v). 



(A.20) 
(A.21) 



(A.22) 
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Using the following identity for the T-functions [91 

(i - v + 1) = £C 

j\{j+P) T(p + v) 



~ r (j-i/ + l) _ T(p)T(u)T(l-u) 



3=0 

it is easy to show that 

U(-v)T(-v) = T(v)U{v) = H(v), (A.24) 
where H{u) is a Hankel matrix: 

(k+j-l)/2 

**M = i + t _ 1 + y - < A - 25 > 

This means that the determinant representation (|A.20p can be rewritten in terms of Toeplitz 
and Hankel matrices: 

f (u,w) = det \l -T^^H^H^T^i-v)] . (A.26) 

Now we have to recall some general results for Toeplitz and Hankel matrices. Let a(z) be a 
smooth function on a unit circle. Its Fourier coefficients are 

2tt 



l]n = 2^/ d9e ' in9a ( eie )- 



We also define the function o such that [a] n = [a]- n . The Hankel and Toeplitz operators 
associated with such a smooth function can be defined as 

Hjk[a] = [a]j+k-u T jk [a] = [a\j_ k , with j,k = 1,2, 

We will use also the following identity for the Hankel and Toeplitz matrices: for any two 

functions smooth on the unit circle, it is easy to demonstrate, using basic properties of the 
Fourier coefficients, that 

T[ab] = T[a]T[b] + H[a]H[b). (A.27) 

Assume that the function a admits the Wiener-Hopf factorization a(z) = a+(z)a-(z) with 
unit constant term, where 

a+(z) = exp |^z fe [lna] fe | , a-(z) = exp j^z" fc [lna]_ fc j . 

Then the corresponding Toeplitz matrix can also be factorized as T[a] = T[a-]T[a + ]. More- 
over, it can be easily shown that, for any function b, the following identities hold: T[ab] = 
T[a-]T[a + b] = T[a-b]T[a+]. 
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The last result we need was obtained by H. Widom for the determinant of products of 
Toeplitz matrices. If a and b are smooth and admit the Wiener-Hopf factorization, the following 
determinant can be computed: 



det [T^la+jTia+b^T- 1 ^]] = exp I ^ k [In a] fc [ln b] 



(A.28) 



.k=l 



Now we can apply these results to compute the determinant (|A.26p . We define the following 
functions: 



n/2 7 n 



W * Z 



n + v 



*(*)= E 



n/2 7 n 



W I Z 



n — v 



(A.29) 



The corresponding Toeplitz and Hankel operators are exactly the ones appearing in ()A.26j) : 

H[a] = -H[b] = H(u), T[b] = T(-u), T[a] = f[y). (A.30) 
It is easy to write an explicit form of these functions. Setting 

w l/2 = ^ we 

obtain 



a(e if> ) 



^ g2iirv ' 



with 9 G [-^,2vr - if)], 



1A qvp iv(Q+2-K-'4>) 

b{e i6 ) = 1 _ e2mv , with G [-27r + V,# 



(A.31) 
(A.32) 



It is important to note that, although these functions are not smooth on the unit circle, all the 
results for the Toeplitz and Hankel matrices mentioned above can nevertheless be used also for 
such functions 



9C|. 



Both functions admit the Wiener-Hopf factorization; it is easy to see that 

[lna] n = In ) <o + ~w n ' 2 (1 - 6 n>0 ) , 

\s1n7rW n 

[\nb] n = In (-^—) ~ ~w n ' 2 (1 - J ft)0 ) • 



Now we can easily compute the determinant (|A.26P : 



det 



(A.33) 



(A.34) 



fo{y,w) = det I + T- 1 [a]H[a]H[b]T- 1 [b] 
where we used (|A.27p . Applying the Wiener-Hopf factorization we rewrite the determinant as 

fo(u,w) = det [r- 1 [a+]T[a + 6_]T- 1 [6_ 
Finally, using (|A.28j) and (|A.33p . we obtain 



OO o 



fo(v, w) = exp < - V — w Tl 
^-^ n 



(!-«;)- 



(A.35) 



n=l 
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A. 3 Summation formula for i ^ 

The general I case can be reduced to the case i = by a simple combinatorial trick which 
we call 'background shift'. It is convenient to consider the sum in (|3.14p as a sum over all 
possible excitations over a Dirac sea (all non-positive integers) which can be constructed as 
either particles (positive integers pj) or holes (non-positive integers 1 — hj) with a fixed difference 
between numbers of particles and holes n p — = i. Equivalently, every term of the sum can 
be considered as an excitation over a shifted Dirac sea (all integers less than or equal to €) with 
particles in the positions £+p k , Pk = 1,2, . . . and holes at l + £ — h k , % = 1,2,... . This second 
parametrization is convenient because there the number of particles is equal to the number of 
holes. 

Before giving the general case we illustrate the background shift in the simplest case £ = 1. 
Let us consider the following sum: 

oo n+1 n 

/im=e e e n^n^ 

n=0 pi< - <p„+i hi <—<h n 3=1 3=1 
p a eN* h a &N* 

n+1 n 

2n n (Pi - Pk) 2 n (hj - h k ) 2 n+1 

sm7w \ j>k j>k tt r z (p k + v) -p-j- T z (h k - v) 

7 n n (Pi + fcfc - 1) 2 fc=i fc=i 1 k) 

3=1 k=l 

We take an arbitrary term of this sum with a given n and given positions of particles and holes. 
We can distinguish two possible situations: p\ = 1 and p\ > 1. In the first case there is no hole 
at point 1 and we define the positions of particles and holes with respect to the shifted Dirac 
sea as follows: 

Pi=Pi+i-l> 3 = l,---,n, hj = hj + l, j = l,...,n. 

Thus we obtain a configuration with n particles and n holes. Substituting this into (|A.36j) . after 
elementary algebra we arrive at 

n+1 n 

n+1 X 2n U(Pj-Pk) 2 U.(hj-h k ) 2 n+1 n 

nw-iTT fe, ( smnu \ j>k j>k yj T^ipk + u) yj T z (h k - v) 
W 11™ n+1 n 11 V2( n A 11 



i=i i=i 



7T / n+1 n 

3=1 k=l 



. 2n 

, SUJ 77// 



i=i i=i 



x det r 




r 2 (p fe + ^ + 1)^(^-1/-!) 

Pj + h k -ij ^ r 2 (p fc )r 2 (/i fc ) 
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In the second case (jp\ > 1) there is a hole at point 1. We set 

Pj=Pj~ 1 i 3 = 1, •••,«+ 1, h\ = l, hj = hj-i + l, j = 2,...,n+l, 

and we obtain a configuration with n + 1 particles and n + 1 holes. Substituting this into (|A.36p . 
we obtain 



n+1 n 
n+1 n , . x 2n I! (Pj ~ Vk? II Oi ~ n+1 



v y n n + % - 1) 2 tm) *-> 1 k> 

j=i k=i 

n+1 n+1 , . 2(n+l) 



x fdew 1 )' Tfr^ + , + i)r^-,-i) 

It is easy to see that, up to a general coefficient T 2 (l + v), we retrieve exactly the same terms 
as in the I = case (jA.ip . It is also easy to observe that the background shift produces a one 
to one correspondence between the configurations of particles and holes in (|A.36p and in ([A.ip . 
up to the replacement of v by v + 1. This leads to a very simple result: 

f 1 (v,w) = r 2 (l + v)f (is + l,w). (A.37) 

We can now use this approach in the case of general I. Let for definiteness £ > 0. We define 
fi{v,w) for < 1 and arbitrary complex v by the following series: 

OO 7l~\-£ TL 

n=0 Pl<-<P„+e hi<-<h n j=l j=l 
p a £N* h a £N* 

Tl~\-£ TL 

. 2n n (Pi - Pfc) 2 n - m 2 n+^ r2 r , , n r2(h , 

smTW \ j>k j>k T j T z (p k + v) -p-r T z (h k - v) 

7T n + i n 11 T -( Pk ) H Mfc) ' 

n n^+^-i) 2 fc=i fc=i 

i=l k=l 

We consider an arbitrary term of this sum with n + £ parameters p k and n parameters hj . We 
introduce a new equivalent set of parameters {p, h} as follows. 

• There is a value m, ^ m ^ £ such that for k ^ m all the parameters p k ^ £, while for 
k > m all the parameters > I. Then for k = m + 1, . . . , n + 1 we define pk-m = Pk — £■ 

• There are £ — m integers ^ ^ ji > j 2 > ■ ■ ■ > je-m ^ 1 such that pk ^ ji, V/c, i. Then we 
define hi = £ — ji + 1, for i = 1, . . . , £ — m. 
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• We define hk+i- m = hk + £, for h = 1, . . . , n. 

As in the case £ = 1 an advantage of this new set of parameters is the fact that there are the 
same numbers (n + £ — m) of parameters {p} and {h}. After some computations (similar to 
those for the I = 1 case but more cumbersome), we obtain 

n+l n 

\2 rr (u . u. \2 



2n II (Pj ~ Pkf II (fy - ^) „+£ r2f \ n r2 



J 1 J 1 1 1. 1 1. >/'./ 

i = l fc=l 



I n+l-m , . s". 



n nipj + hk-i) 2 k=i y ' k=i 

n+£-m , . N 2(n+l-m) 

/ cin tj; \ 

3=1 3=1 

^ + m Pi + ht-i) £ P"(p*)r a (fc*) 

It is easy to note that there is a one to one correspondence between all the possible configurations 
{p, h} and {p, /i}. Hence we retrieve the sum (|A.lj) in terms of the parameters p and /i with a 
simple prefactor: 



\3=1 



" 1)/2 ( n r V + j) ] (l-w)-^ . (A.39) 
\j=i 



For negative £ the computation can be performed in a similar way leading to the following 
result: 

f e (u,w) = w^ 2 y}r-> + l-j)j (l-^)-^) 2 . (A.40) 

The results for positive and negative values of £ can be written together using the Barnes 
G-function: 

Mu,v,) = ^^ ^(t+y (1 " W) ' {V+t? - (A - 41) 

Now setting w = e 2mx / L ~ £ ; with < x < L, and taking the limit e — > +0, we complete the 
proof of equation ()3.14[) . 
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